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Abstract. We present the results of magnetic properties and entanglement of the distorted diamond chain 
model for azurite using pure quantum exchange interactions. The magnetic properties and concurrence 
as a measure of pairwise thermal entanglement have been studied by means of variational mean-field like 
treatment based on Gibbs-Bogoliubov inequality. Such a system can be considered as an approximation of 
the natural material azurite, Cu3(C03)2(OH)2. For values of exchange parameters, which are taken from 
experimental results, we study the thermodynamic properties, such as azurite specific heat and magnetic 
susceptibility. We also have studied the thermal entanglement properties and magnetization plateau of the 
distorted diamond chain model for azurite. 



1 Introduction 

Copper oxide materials as low-dimensional magnetic sys- 
tems are interesting subjects to investigate because of the 
new physics that can arise in low temperatures. The nat- 

\ ural material azurite Cu3(C03)2(OH)2 is the one of such 
copper oxide material which has been the subject of de- 
bates in recent years |1I2I3I4I6I7I8| . 

Moreover, azurite can be considered as one of the first 

" experimental realizations of the ID distorted {J\ ^ J3) 
diamond chain model (see Figure [TJ for the detailed struc- 

■ ture of azurite see for example [I]). It shows antiferromag- 
nctic behavior at temperatures below Necl temperature 
1.9 K [S]. Both the experimental analysis and theoreti- 
cal modeling have been performed, but there is no clarity 
about exact values of the relative exchange interactions 
between azurite's copper ions [1I2I3I4I6I7I8J . On the other 
hand, there is significant evidence that the magnitude of 
the exchange interaction J2 is the greatest one, forming 
dimcrs and monomers of Cu- ions |2I4I7) . In [5] an effective 
generalized spin- 1/2 diamond chain model has been sug- 
gested with a dominant antiferromagnetic Cu-ions of the 
dimer coupling J2, two antiferromagnetic Cu-ions of the 
dimer- monomer couplings J\ and J3, as well as a signif- 
icant direct monomer- monomer of Cu-ions J m exchange 
which explains a broad range of experiments on azurite 
and resolves the existing controversies. 

The most noticeable feature is that there is a large 
magnetization plateau at 1/3 of the saturation magne- 
tization [3], which extends from 11 to 30 T, when the 
magnetic field is applied in the perpendicular to the chain 
axis. Such 1/3 plateau is usually associated with classi- 
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Fig. 1. Distorted diamond chain. 



cal up-up-down (uud) type of spin arrangement (or with 
a quantum state which has a classical analogue), while 
the 1/3 plateau of azurite is proposed to be of fundamen- 
tally different, OOu type, where the dominant J2 coupling 
ensures that the two "dimer" spins on the Cu2 sites (see 
Figure [I]) are in a singlet state, while the third "monomer" 
(Cui) spin is completely polarized by the magnetic field. 
This state is based on the presence of a singlet and is a 
pure quantum nature without a classical analogue. There- 
fore azurite is the first successful candidate for describing 
a 1/3 quantum plateau state [10111) . Another characteris- 
tic feature is that there is almost a direct transition from 
the plateau of saturation, which may be an interpretation 
as a remnant of the so-called localized magnon-members 
present in a perfect chain of diamond J12I13] . 

The phenomenon of magnetization plateau has been 
studied during the past decade both experimentally and 
theoretically The plateau may exist in the magnetization 
curves of quantum spin systems in the case of a strong 
magnetic external field at low temperatures. Magnetiza- 
tion plateaus appear in a wide range of models on chains, 
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ladders, hierarchical lattices and theoretically analysed 
by dynamical, transfer matrix approaches as well as by 
the exact diagonalization in clusters [14115116117118111?] . 
In this paper, we obtain the magnetization plateau in ID 
diamond chain using variational mean-field like treatment, 
based on Gibbs-B ogoliub ov inequality J20I21I22I23I24I25I26] . 

Entanglement [27I30J has gained renewed interest with 
the development of quantum information science. Entan- 
gled states constitute a valuable resource in quantum in- 
formation processing |31l32j for example the predicted ca- 
pabilities of quantum computers rely on entanglement. 
Numerous different methods of entanglement measuring 
have been proposed for its quantification [27] . In this pa- 
per we use concurrence |33I34] as entanglement measure of 
the spin- 1/2 system. Geometrical frustration and thermal 
entanglement (concurrence) of spin- 1/2 Ising-Heisenberg 
model on a symmetrical diamond chain was studied in 
[35] . In this paper we study the concurrence of spin- 1/2 
Heisenberg model on a distorted diamond chain as the 
approximation of natural mineral azurite. 

This paper is organized in the following way. In sec- 
tion 2 the variational mean-field formalism based on the 
Gibbs-Bogoliubov inequality is applied to the Heisenberg 
model on the distorted diamond chain. In section 3 we in- 
vestigate the magnetic properties of the system and com- 
pare the obtained results with the experimental data of the 
magnetization, specific heat and magnetic susceptibility. 
In section 4 the concurrence of the system is calculated. 
Some conclusive remarks are given in the last section. 




Fig. 2. Each Ai cluster consists of one rectangle of S\, 
S 2 , S 3 , S\ sites (grey rectangle) and dimer of S 2 ,, S\, 
sites (bold line). 



where H is the real Hamiltonian which describes the sys- 
tem and Hq is the trial one. F and Fq are free energies 
corresponding to H and Hq respectively and (. . .) de- 
notes the thermal average over the ensemble defined by 

By introducing trial Hamiltonian for our model con- 
taining unknown variational parameters one can minimize 
right hand side of Gibbs-Bogoliubov inequality (J3J and get 
the values of those parameters. 

We introduce a trial Hamiltonian Hq as a set of nonin- 
teracting clusters of two types in the external self-consistent 
field: the rectangle and the line (see Figure [2]) 



Hq = J^(Hq^+Hq(% 



(5) 



2 Gibbs-Bogoliubov approach 

We use a spin-1/2 Heisenberg model. The Hamiltonian of 
the Heisenberg model is 



A, 

where the indexes of summation Ai 
tcracting clusters and 



label different nonin- 



H {i) 

n 



= \'{s 2 ,si,)- l2 sy- 



- 73SI 2 



i±s\ 
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H — 2_^ JijSiSj — y gusBSj, 
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(\\ where A, A' and 7j,7^ are the variational parameters. 

The real Hamiltonian ([2]) can be represented in the 



following form: 



# = £#« 



where Si are the spin-1/2 operators, Jij is the exchange 

interaction constants connecting sites i and j, B is the 

value of the external magnetic field, g - the gyromagnetic 

ratio and fi B - the Bohr magneton. The Hamiltonian for where H (i) hag thg following f on 

the distorted diamond chain can be written as 



(7) 
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where h = g/j,BB z and the g is set to 2,06 [55] and 



on = S2S3 + 5 m S\S\ + S2(S\S % 3 + S l 2 S\) 
+ S 3 (S\S 2 + SIS\). 



(3) 



•h (S2/S3/J - 

- h(sy + sy) + ^{sisi" 1 + s\s\ +1 ) + 

+ Ji(SlS*7 x + S\S 2 ,) + J 3 (S1S' 2 7 1 + S\Sl,), (8) 

Gibbs-Bogoliubov inequality Q can be rewritten now as 
follows: 

/ W < / (i) + (J2 - A) (Oi) + (J2 - A') <S 2 ,S 3 ,) + 



(9) 



where S 2 = J\j J 2l 63 = J3/J2 and S m = J m /J 2 - Here and 
further exchange parameters (Ji, J 2 ,J 3 ) and the magnetic 
field h are taken in Boltzmann constant scaling i.e. Boltz- 
mann constant is set to be fcg = 1. 

Gibbs-Bogoliubov inequality [20I21I22I23I24I25I26J states 
that for free energy (Helmholtz potential) of the system 

we have where f( % > = F/N and Jq = Fq/N are free energies of 

F < Fq + (H — Hq) q , (4) the one cluster (N is number of clusters) and we denote by 



Ji (mim 3 + m^m 2 ) + J 3 (m\m 2 + 77147713)+ 

4 3 

J m (mim 4 ) - J2 (h~ 1]) m j - J2 ( h ~ l'j) m 'i 

J'=l J'=2 
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Fig. 3. Magnetization m of the cluster versus external Fig. 4. Magnetizations m 1 ,m 2l m 3 and m 4 versus ex- 
magnetic field h (Tcsla) for J 2 = 33A", J x = 15.51 K, tcrnal magnetic field h (Tcsla) for J 2 = 33A', J x = 
J 3 = 6.93 A, J m = 4.62 A at and T=1.3 K. 15.51 K, J 3 = 6.93 A, J m = 4.6 K. at T=0.8 K. 



nij = (Sj Z ) n and m'-, = \Sji z ) n the magnetisations of A* 
cluster. Here we take into account that spins of different 
clusters are statistically independent, i.e. (SjSff ) = 

mjTn' k . For given h, Ji one must perform the minimization 
of r.h.s and obtain the values for 7, A and 7', A'. Minimiz- 
ing the right hand side of © in order to 7., , A and 7J , A' 



and using the fact, that 



dX 



= (<Xi) and 



we obtain the following values for the variational param- 
eters: 

A = A' = J 2 , 
72 = 73 = h, 

71 = h - Jim' 3 - J 3 m' 2 - J m rrii, 
74 = h - J 3 m' 3 - Jim' 2 - J m mi, 

7 2 = h — J1TO4 — J3ITI1, 

73 = h — J3?7i4 — J\m\. (10) 

Using this values and the trial Hamiltonian one can cal- 
culate the value of any thermodynamical function of our 
system for the fixed h, Ji- 



3 Magnetisation, specific heat and 
susceptibility 

The results of the previous section can be used for in- 
vestigation of the magnetic properties of the model. The 
magnetization of arbitrary site m v of cluster Aj is defined 
as 

_ Tr(S i v 'e-< > / T ) 



(11) 



(i) 



where S % v z is corresponding spin operator, H is the Hamil- 
tonian ([6]) and Z is the corresponding partition function 
of the cluster. To obtain all six magnetizations (m 2 ,m 3 , 
mi, . . . 7714) one must insert the values of the variation pa- 



rameters pOj) into the equation (fTTj) and solve the result- 
ing system of equations for the fixed h, Ji. 



mi = fi(m 2 ,m' 3 ,mi,mi) 



mi 



/4K: 
fsimi; 



m 3 , mi, 7714J 

m 4 ) 

m 4 ) 



(12) 



The calculations of our paper is based on the effective 
diamond chain with values of coupling constants taken 
from [2]: 

J 2 = 33A; 5 2 = 15.51/33; S 3 = 6.93/33; S m = 4.62/33. 

(13) 
The dependence of the average magnetization (for clus- 



ter) 



id 



mi + m 2 + m 3 + m^ + m' 2 + m 3 
6 



(14) 



from external magnetic field, calculated from {T2J , is shown 
in Figure |3l As it can be seen from the Figure [3] the 1/3 
magnetisation plateau at T = 1.3K extends from 11 T 
to 29 T interval, while the experimental data |2I4| gave 
11 T - 30 T interval. 

As it mentioned above, azurite is a good candidate 
to exhibit OOu type 1/3 plateau state jlOlllj . The plots 
in Figure 2] illustrate magnetisation versus magnetic field 
dependencies for different sites obtained using mean-field 
approach. As it can be seen from the Figure [4] dimers are 
essentially in the singlet state (dashed lines) whereas the 
single "monomer" spins are almost fully polarized in the 
1/3 plateau (solid lines). The observed polarization is in- 
compatible with a uud type of plateau, in which the dimer 
spins are strongly polarized. We find that dimer spins are 
about 2.5% polarized while the numerical calculation for 
the ideal diamond chain Heisenberg model gives 2.7% po- 
larisation [5], while experiments give 10% [TU] . 

The magnetic susceptibility is defined as follows: 



Xo 



dm 
~dh 



(15) 



h=0 
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Fig. 5. Zero field magnetic susceptibility versus temper- 
ature. 




10 20 30 T 

Fig. 6. Zero field specific heat versus temperature. 



The magnetic susceptibility measurements of azurite per- 
formed in j3| and it was found double-pcak-like structure 
in the magnetic susceptibility curve, namely, in the tem- 
perature dependence of the magnetic susceptibilities the 
sharp peak appears at 5K and the rounded peak is ob- 
served at 23K . The Figure [5] shows the temperature de- 
pendence of magnetic susceptibility obtained by our cal- 
culations using (fT2")) and (fT^|) . It has first sharp peak at 
AAK. 

Analogous to the magnetic susceptibility, we also have 
calculated specific heat: 



O 



T fd 2 fo 
\dT 2 



(16) 



h=0 



The specific heat measurements for azurite are performed 
in [3] and a sharp peak is observed at Tn = 1.8K signal- 
ing an ordering transition and two anomalies have been 
observed in the specific heat at T = 4K and T = 18K 
(415) . The first peak is out of reach of a one-dimensional 
model. Our calculations gave a low-temperature peak for 
H = at T = 3K and the second peak at 12 K (Fig- 
ure [6]). The obtained double-peak-like structure of specific 
heat reproduce the important features of the experimen- 
tal results 03]. Density functional theory (DFT) also 
gives the double-pcak-like structure for the specific heat. 
In this paper we have reproduced the important features 
of the experimental results in the specific heat properties 
of azurite theoretically. 



4 Concurrence and thermal entanglement 

The mean-field like treatment transforms many- body sys- 
tem to the set of clusters in the effective self-consistent 
magnetic field. This allows to study, in particular, ther- 
mal entanglement properties of the system. The entangle- 
ment of formation [29130] was proposed to quantify the 
entanglement of a bipartite system. Unfortunately the en- 
tanglement of formation extremely difficult to calculate, 
in general. However, in the special case of two spin 1/2 
particles an analytical expression |33I34| can be obtained 



for the entanglement of formation of any density matrix 
of two spin 1/2 particles: 



£F = „ ( <i±4ZE) ) , 



(17) 



where 



H(x) = -x\og 2 (x) -(1-x) log 2 (l - x), (18) 



and C is the quantity called concurrence |33l34j . In the 
case of diamond chain the concurrence Cij corresponding 
to the pair (i,j) has the following form: 



a 



i -j 



;{Ai - X-2 - A 3 - A 4 ,0}, 



(19) 



where A& are the square roots of the eigenvalues (Ai is the 
maximal one) of the matrix 

^fti'W^D'PrW^J), (20) 

where pij is the reduced density matrix for (i,j) pair. For 
rectangle cluster the density matrix is 



P = 



1 16 



'l^»)(^»|. 



(21) 



where Z is the partition function of the rectangle and ipi 
and Ei are eigenvectors and eigenvalues of the Hamilto- 

nian Hq J6J respectively. The reduced density matrix of 
(i,j) pair pij can be calculated as 

Pij = z2(4 l \p\$)> ( 22 ) 



where </>?f is the p-th basis vector for (k,l) remaining pair 
of sites of rectangle. After the calculations it has the fol- 
lowing form 

'u 0\ 

v J 
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Fig. 7. The concurrences Cj 2 ,Cj i: C,j 3 and Cj m versus i q 
external magnetic field h (Tcsla) at T = 0.1K for J 2 = 
33AT, Ji = 15.51 K, J 3 = 6.93 K, J m = 4.62 X. 



where u, w, y and v are some functions of variables 7, A, 7', A' 
and T. Using flU]), (gDJ and &$ one can find the following 
expression for the concurrence Cij 




Cij = 2max{\y\ — ^/uv, 0}. 



(24) 



In Figurcjjjis shown the dependence of concurrences 6*2,3 = 
W 2 , 61,2 = C3.4 — Cj 3 , 6*1,3 = 62,4 = Cj x and Ci,4 = 
Cj m on the magnetic field. The behavior of the concur- 
rence can be used to analyze spin phases of azurite. As 
it can be seen from the Figure [7J there is three regions 
in magnetic field axes with different ground states. For 
lower value of external magnetic field the opposite spins 
(5*2 and S3) in diamond cluster is highly entangled. The 
neighboring spins with lower coupling constant are not 
entangled. For higher values of h when the magnetization 
has a plateau the entanglement of (Si, £4) pair is almost 
zero, i. e. practically unentangled, while the (52, 5*3) pair 
is almost fully entangled. The concurrence of the neigh- 
boring spins on the plateau is small, comparing to (5*2, 6*3) 
pair, moreover the neighboring spins with lower coupling 
constant ( J3) are unentangled. 

Now, we consider the dependence of the concurrence 
on temperature. Figure [5] (a) shows the temperature de- 
pendence of the concurrence for different pairs of diamond 
chain at small value of the magnetic field (h = IT). The 
neighboring spins with lower coupling constant (J3) stay 
unentangled with increasing temperature while the con- 
currence of the bigger one (Ji) decreases with tempera- 
ture to 4.5-ftT, where the entanglement vanishes. The con- 
currences of J m pair J\ behave similar and vanish almost 
at the same temperature, while the dominant Ji pair stays 
entangled at higher temperatures (until 28.&T). Almost the 
same behaviour shows the temperature dependence of the 
concurrence for different pairs of diamond chain at plateau 
phase h = 18T (Figure[5](b)). The concurrences for J m , J 3 
and J\ decrease with temperature and vanish sequential 
between AK and IK while the concurrence for J2 pair 
stays entangled for the higher temperatures and vanishes 
at the same temperature as for small values of magnetic 
field (2SK). 
As it can be seen from the Figure [JJ and [5)3 in the plateau 




Fig. 8. The concurrences Cj 2 , Cj 1 , Cj 3 and Cj m versus 
temperature for a) h = 0T and b) h = 18T. 



state dimers are almost fully entangled (lines labeled by 
Cj 2 in the figures) whereas the monomer spins are weakly 
entangled (lines labeled by Cj 1 and Cj 3 in the figures). 
The observed entanglements is incompatible with a (uud) 
type of plateau, and confirm the proposed (OOw) nature of 
the plateau. 

Now we revert to the Figure [5] (a) to notice that com- 
parison of the Figure [S] (a) with the Figure H] shows that 
the C(J2) has a peak at nearly T = 5K and it is located 
between two peaks of the specific heat. Roughly such a 
behavior can be understood as follows. 
As a result of interaction between the horizontal (J m ) and 
vertical ( J2) dimers and also as a result of an asymmetry 
(J\ > J3), decreasing of the concurrences (in comparison 
with non-interacting case) of these dimers at zero temper- 
ature is observed. As temperature increases, energy "ac- 
cumulates" in the horizontal dimer at first and in the ver- 
tical dimer later, which causes the double peak structure 
in the specific heat picture (see Figure [5]). During the pro- 
cess of " energy accumulation" in the sites of the horizontal 
dimer (the first peak region) destruction of a quantum cor- 
relations between them takes place as a result of thermal 
fluctuations (so the concurrence C(J m ) is decreasing, see 
Figure [8] (a)). And also destruction of a quantum corre- 
lations between the sites of horizontal dimer and vertical 
one occurs. As a consequence one can see an increasing 
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of C(J2) from T = OK to T = 5K (without above men- 
tioned couplings (non-interacting dimers case) C( J2) gets 
its maximal value equals to 1 in this region) . Further tem- 
perature increasing brings destruction of the quantum cor- 
relations between the sites of the vertical dimer, which is 
the result of "energy accumulation" on these sites (and 
as a consequence - decreasing of 0(h) and increasing the 
specific heat to the second peak). 

Varying values J\, J2, J3, J m of the coupling constants, 
brings to analogical picture, except the symmetric case 
where J\ = J3. In the symmetric case the concurrences 
C(Ji) = C(J m ) = 1 (at zero temperature) and are de- 
creasing with the temperature and as fast as are higher 
the values of the J\ = J3 coupling constants. Closer is the 
diamond to the symmetric case, worse is the appearance 
of the peak structure in the C( J2) picture. 



5 Conclusions 

In this paper using Hcisenberg model the distorted dia- 
mond chain was studied as approximation for natural ma- 
terial, azurite Cu 3 (C03)2(OH) 2 . The magnetic properties 
and concurrence as a measure of pairwise thermal entan- 
glement of the system was studied by means of variational 
mean-field like treatment based on Gibbs-Bogoliubov in- 
equality. In our approach for the values of exchange pa- 
rameters taken from theoretical and experimental results 
we have obtained the 1/3 magnetization plateau with cor- 
rect critical values of magnetic field, moreover this plateau 
is caused by OOw type plateau state. We also studied the 
thermal entanglement properties of the distorted diamond 
chain and drew a parallel between them and the specific 
heat ones. 
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